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ABSTRACT 

This report summarises work done on "Vehicle Sys- 
tems and Payload Requirements Evaluation" during 
the five years of the study. Reference is made to 
specific previous Interim Scientific Reports. Work 
done since the last Interim Scientific Report is 
reported on in more detail. This includes the 
analysis of Redundant Systems. An investigation of 
potential problems with the Gravitational Red Shift 
mission, and development of the Interactive Graphic 
Orbit Selection computer program. 
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VEHICLE SYSTEMS AND PAYLOAD REQUIREMENTS EVALUATION 
(October 1970 - December 1975) 

by 

F. G, Rea, J. L. Plttenger, R. J« Conlon, and J, D« Allen 


INTRODUCTION 


This report presents the results of the work on "Vehicle Systems 
and Payload Requirements Evaluation" (VSPRE) for NASA Marshall Space Flight 
Center (MSFC) under Contract Number NAS 8-26491 between October, 1970 and 
December, 1975. Background information and a summary of the study results 
are presented for the entire project along with detailed technical discussion, 
recommendations, and conclusions of the past 18 months. Details of the first 
four years have been published earlier in Interim Reports (Reference 1, 2, and 
3 *). 


Program Oblectives 

The objective of this program was to develop techniques for 
identifying launch vehicle system requirements for NASA automated space 
missions, and to conduct analyses to assist NASA OSS in the managemeat of its 
related Supporting Research and Technology (SR&T) programs. AchieveTiicnt of 
this overall objective required development of various ca puter programs. 
These include extensions of previously developed system analysis and evalua- 
tion programs for the analysis of a oroader scope of vehicle systems. In 
addition, these programs were applied to the analysis of several NASA mission 
astrionics requirements. For these missions, trade-off evaluation of com- 
petitive systems and analysis of operational requirements were performed. 

Project Tasks 


The efforts of ttiis study have been conducted under specific 
tasks. The contract modification of April, 1974 redefined the task numbers 
from those in use since 1970. For convenience in referring to earlier re- 
ports, Table 1 shows the old and new task number definitions. All tasks 
discussed in this report are under the new numbers. 


^References are listed at the end of this report. 
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TABLE 1. VSPRE PROJECT TASKS 


Task Description 

Old Task No. 
(Prior to April, 1974) 

New Task No. 
(After April, 1974) 

Update and Maintain 

1 

1 

Computer Codes 

Investigate KEAO 

2 


Astrionlcs Requirements 

Investigate Astrionlcs 

3 

2 

for OSS Missions 

Investigate Scout 

4 

3 

Astrionlcs Requirements 

Perfoni: Mission 


4 

Requirejpents Assessments 



** This effort was begun in 1973 

under Old Task 3* 
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SUMMARY 


The following paragraphs summarize the effort on each of the 
project tasks. Detailed discussions appear in the Technical Discussion 
Section of thlc report. 

Update and Maintain Computer Codes (Ta s k 1) 

Several computer programs have been developed under this task 
and delivered to MSFC. Most significant of these are; 

• EOMP-I (Earth Orbit Mission Program - I) which performs 
linear error analysis of launch vehicle dispersions for 
both vehicle and navigation system factors, and 

• IGOS (Interactive Graphic Orbit Selection) program which 
allows the user to select orbits which satisfy mission 
requirements and to evaluate the necessary injection accuracy. 
This program was written as part of the Task 4 investigation 
of mission requirements for astrlonlcs improvements. 


HEAP 

The High Energy AstronOTiical Observatory (HEAP) mission was 
investigated undar the original VSPRE Task 2. At that time it was proposed 
that HEAO be launched on a Titan HID with a Lockheed Orbit Adjust Stage (OAS) 
for final orbit placement. VSPRE efforts were; 

• Calculation of the Titan illD open loop guidance system 
injection errors, 

• Analysis of the probability of obtaining sufficient 
parking orbit life time to adequately track the 
vehicle, and 

• Development of an optional OAS thrusting policy. 

Results of these studies (References 1, 2 and 3) indicated the Titan IITD/OAS 
vehicle was feasible. 

Subsequently the HEAO payload weight was reduced, and the payload 

was reassigned to an Atlas vehicle. 
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Inv€»tlaate Astrlonlcs for OSS Mlaalong (Task 2) 

A number of specific analysis studies have been directed by 
the COR (Contracting Officer's Representative). Astrlonlcs requirements 
for HEAO and LAGEOS, have been discussed in previous repot t» (Reference 1, 

2 and 3). The Red Shift mission requirements for Scout are discussed 
under Task 3 in this report. 

Recent Task 2 activity has concentrated on an evaluation of the 
requirements for a high reliability computer, specifically the ARMMS 
(Automatically Reconfigurable Modular Multiprocessor Sy .em) c«vnputer 
under development at MSFC. The evaluation of ARMMS resulted in the writing 
of a preliminary version of a generalized computer program, ABBACUS for the 
analysis of a reconfigurable astrlonlcs system. 

Investigate Scout Astrlonlcs Reauirements (Task 3) 

Several potential modifications to Scout astrlonlcs have been 
evaluated and the results previously published in References 2 and 3. 

No conclusive recommendations regarding the need for the modifications could 
be made due to a lack of firm mission requirements for increased Injection 
accuracy. 


Most recently. Task 3 effort concentrated on the application of 
Scout for the Gravitational Red Shift mission. This mission is a sub-orbital 
launch with the payload attached to the spinning, depleted ^ fourth stage. 

The results of the study Indicate no problems are encountered due to the 
operation of the payload while attached to the spinning fourth stage. 

Perform Mission Requirements Assignments (Task 4) 

When establishing mission requirements, the planners frequently 
stated, as a requirement, the known accuracy of the vehicle they assume 
would be assigned to their mission. In other cases, arbitrarily small 
error*- are specified to minimize the chance of failure without regard to 
the Implied Increased astrlonlcs costs. 

This situation led to the development of the Interactive Graphics 
Orbit Selection (ICOS) programs. IGOS displays both the mission requirements 
(orbits which met the requirements) and the various launch vehicle capabilities 
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(payload weight and Injection accuracy). The program serves as the basis 
for concise dialog between mission and vehicle oriented specialists. IGOS 
is a tc''l for either orbit selection (mission design)* or for the evalua- 
tion of launch vehicle accuracy requirements. In addition to the details 
presented in the Technical Discussion, a users manual and a sample work 
session are presented in Appendices A and 6. 


BACKGROUND 
Related Prior Work 


The astrlonlcs/avlonics evaluation techniques extended under 
the contract were originally developed by Battelle for NASA/OSSA under 
contract with the NASA Electronics Research Center (ERC), Contract Number 
NAS12-550, and for the U. S. Air Force (IBAF) Avionics Laboratory, Contract 
Number F33615-69-C-1402. The ERC contract was originally Initiated to develop 
an evaluation technique for inertial sensors to be included in strapdown 
inertial navigation systems. To effectively evaluate sensors for strapdown 
application, it was determined that a more comprehensive systems level study 
including other astrionics subsystems was necessary. 

The systems evaluation techniques developed under the ERC contract 
were extended to cover avionics under the USAF Avionics Laboratory contract. 

In analyzing Integrated avionics systems, the assumption that each item of 
hardware is associated with a particular function (such as navigation or target 
acquisition) is too restrictive. For this reason, a more general technique 
was developed in which a total system is considered as an assembly of sub- 
systems. The operational objectives of the combined system are then defined 
in terms of functions such as navigation and flight control. Any hardware 
subsystem might be involved in satisfying one or more functions. The ERC 
and Avionics Laboratory studies are described in more detail in the following 
paragraphs . 



Techniques for Evaluatln 2 Astfionlca 
for Interplanetary Missions 


Techniques for evaluating astrlonlcs systems for interplanetary 
probe missions were developed under Contract NAS 12-550 over a 3 year period. 
The objective of the study was to develop a technique for evaluating strapdown 
guidance systems. The technique utilizes system parameters describing the 
reliability, power, weight, and performance of the strapdown guidance system 
in determining a measure (or index) of guidance system performance. This 
work is discussed in detail in References 4 through 10. 

During the work performed for ERC it was assumed that the astrlonlcs 
being considered were an Integral part of the spacecraft as shown in Figure 1. 
In using this program, it would not be necessary for the astrlonlcs to be an 
integral part of the spacecraft if only launch vehicle navigation and guidance 
were being considered. In this case, certain subsystems such as the space- 
craft attitude control would not be considered in the analyses. 

Techniques for evaluating astrlonlcs for interplanetary missions 
require definition of a baseline set of Inertial sensors, computers, and 
other onboard hardware items that are not designed internally by the program. 

In addition, parameters that describe the mission must be provided to generate 
the appropriate flight trajectory. 

One limitation of these evaluation techniques is that no provision 
is made for consideration of backup modes for the various functions, e.<:., 
navigation. The user must specify those subsystems he wishes to use to perform 
each function. The evaluation yields measures of the penalty for Including 
those subsystems. Work initiated under USAF contract in 1969 allows considera- 
tion of backup modes for various functions. This work is described in the 
following section. 

Techniques for Evaluating 
Integrated Avionics 

The objective of this study (performed under Contract Number 
F33615-C-1402) was to develop techniques and associated computer programs 
for evaluating the effectiveness of integrated avionics for USAF aircrait 
on appropriate missions. Further discussions of this work are contained in 
References 11 and 12. 




* Deep Space Instrumentation Facility 


FIGURE 1. BLOCK DIAGRAM OF INTEGRATED ASTRIONICS 











Battelle's Navigation Error 
Analysis Program (NEAP) 


The technique of performing linear sensitivity analysis used 
under the ERC contract to determine strapdown navigation errors was genera- 
lized to include all possible navigation mechanizations and extended to 
cover open-loop vehicle dispersions under a Battelle-funded study described 
in detail in Reference 13. Experience gained from this Battelle-funded 
study was used in writing the Earth Orbital Mission Program-I (EOMP-I). 

Concurrent Related Projects 

Concurrently with performance of Contract Number NAS 8-26491, 
Battelle personnel were Involved in two related projects funded by NASA 
OSS which provide valuable data and practical experience. 

NASA Launch Vehicle Project (NLVP) 

This project (Contract Number NASw-2018) provides extensive launch- 
vehicle and trajectory data for present and future NASA expendable launch 
vehicles. Model formulations and computer programs are structured to maintain 
compatibility with NLVP data. In addition, consultations with NLVP personnel 
on the subjects of propulsion and vehicle capability provided valuable support 
for this study. 

Delta Inertial Guidance 
System (DIGS) Project 

Battelle assisted the Delta Project Office (DPO, Goddard Space Flight 
Center) in the technical evaluation of the Delta inertial guidance system. 

Work under this contract (NAS 7-786) has Included the development of a FORTRAN 
program to simulate the Delta launch vehicle, its control systems, and on-board 
computer software. This provided VSPRE personnel with a very detailed des- 
cription of the NASA OSS launch vehicle, onboard systems, and software. As 
discussed in Reference 3, the DIGS project has also provided an opportunity for 
an independent check of results obtained by using EOMF-I. 
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TECHNICAL DISCUSSICTt 


Tt^ie Technical Discussion section of this report Is divided 
Into discussions of the four current project tasks. The recotmnendations 
resulting ;rom each of the tasks are presented separately In the 
Reconr.teiid>;cions section later In this report. 

Tafik \ Development and Maintenance of Technlouea and Computer Cod es 
for the Evaluation and Analysis of Astrlonlca Systems 

Several computer codes are used In the analysis of astrlonlcs 
systems, "hese programs are summarized In Table 2. The table Includes 
reference; to prior reports for programs developed and/or used during 
earlier reporting periods. 

The Interactive Graphics Orbit Selection (IGOS) program has 
been written during this reporting period and has recently been delivered 
to MSEC for operation on the UNIVAC 1108. This program operates inter- 
actively from a remote graphics terminal such as a Tektronix 4010/4012. 

A complete users manual is included in Appendix A and a sample 
work session is included in Appendix B. Technical details of the program 
are des''ribed under Task 4 of this report. 



I § 

10 

I ■’ 

; . i 

i 

' i 

* * 

1ABLE 2. VSPRE COMPUTER CODES 


1 

f 

PROGR.^ NAME 

DESCRIPTION 

STATUS* 


EOMP-I 

Generates Launch Vehicle injection disper*^ 
slons (References 1, and 2) 

Delivered to MSEC. 
Converted to 1108 


EOMP-II 

Compares Launch Vehicle Injection disper- 
sions to mission requirements (Reference 2) 

Inactive due to ^ 

lack of adequate 
mission requirements. 

"A 

MONTE CARLO 

Computes statisfled distribution of 
orbital elements from Guassian injection 
covariance (References 2, and 3) 

Delivered to MSEC. 
Operational on 
1108 


Performance Code 

Generates optimal nominal trajectory from 
vehicle weight statement (Reference 3) 

Not developed on 
VSPRE contract. 

i. 

-h 

ARBACUS 

Modular Redundant System Reliability 
Analysis (Task 2 of this report) 

Preliminary version 
only. Development 
dropped. 


SPUN STAGE 
SIMULATION 

Simulates a generalized spinning body with 
mass flow and thrusting. (Reference 3 
and Task 3 of this report) 

Preliminary program 
for specific applies-- 
tions. Not delivered. 


IGOS 

Interactive Graphics Orbit Selection 
(IGOS) program. Generates mission 
accuracy requirements. (Task 4 of this 
report) 

Delivered to MSFC. 
Operational on 1108 


* All programs are operational on Battelle's CDC 6400/Cyber computer system. 
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Task 2 Investigate Aatrlonlcs Requirements for OSS Mtaslone 

High reliability modular computer development has received considerable 
funding under the NASA OSS Launch Vehicle and Propulsion Programs SR&T program. 

In particular, the Automatically Reconflgurable Modular Multiprocessor Systems 
(ARMMS) project at NSFC was the prime reason for most of the effort. The VSPRE 
project team was directed by the COR to review the ARMMS project and to evaluate 
the concept particularly In terms of Its impact on vehicle reliability. VSPRE 
memoranda (References 14 and IS) were published as the result of attendance at 
an ARMMS design review and a study of ARM^C! documents. 

Modular Redundant System Reliability 


The evaluation of modular redundant systems is a more general subject than 
the ARMMS reliability study. The Insight gained and the analytic tools developed for 
the ARMMS-related Investigations are applicable to many launch vehicle systems and 
subsystems. Including redundant astrlonlcs hardware. Therefore, the analysis was 
carried through to the development of a computer code. Its use and its integration 
with other techniques are presented below. 


Development of an Analytic Method 

In order to evaluate the reliability of a modular redundant system, it Is 
necessary to know the loading on the various pieces of hardware. A nominal mission 
is deterministic. If no hardware failure.? occur, then the on-off history of every 
component Is known a priori. Thus, the on and off failure rates can be used as 
appropriate In a system simulation. Analytical problems arise because. In any real 
situation, the on-off histories are not deterministic. Generally, the evaluation 
of reliability must be accomplished in a multi-dimensional space. The on-off 
histories are correlated In a very complex 'ay. Every component state Is 
propagated as a function of the dynamic states of all the system modules. This 
kind of problem can be handled by means of convolution Integrals In which the final 
probabilities of failure for the spare modules are dependent upon both the on and 
off failure rates of the modules vrtilch are nominally planned for use. An active 
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module which has a very high failure rate would tend to decrease the reliability 
of a back-up module near the end of the mission, since there Is a higher probability 
that It has been used at that point. 

The predicted reliability of any module In the system Is dependent upon 
the previous history and records of failures and use rates of every module In the 
system. The analysis necessary to handle this problem la exceedingly complex, and 
one Is finally forced to use a Monte Carlo type of approach. For many appllca* 
tlons this would be too expensive and, therefore, some simplifying assumptions have 
to be made. 

During the course of the ARMMS study, one of the parameters of interest 
was the ratio between the active and dormant failure rates of modules. If, as has 
been indicated by a recent study (Reference 16), failure rates for active and dormant 
modules differ not by orders of magnitude but only by a factor of two or less, then - 
In terms of numbers one normally deals with and accuracies which are required for 
reliability computations - it is not necessary to include in the model separate 
failure rates for the modules It. their active and passive states. The Inclusion 
of an extra module might Improve reliability by two or three orders of magnitude; 
whereas, accounting for different on-off failure rates would not affect the calcu- 
lation nearly this much. 

Another problem common to this area of reliability analysis is the large 

number of individual component reliability states which must be propagated. For a 

N 

system with N modules, there are 2 different failure states the system may be in. 

For a system of only 20 modules, this number is already in the neighborhood of a 
million, and certainly far beyond anything that could reasonably be treated on an 
independent state basis. It is necessary to view the problem in such 5 way that 
fewe- states are encountered. This is possible because, for purposes of the 
re'tabil'ty analysis, we only need to know how many of a given type of module 
have fal’-‘d - not specifically which ones. It Is assumed that the executive system 
will correctly keep track of which modules are operable. Let us define, for each 
set of N modules of a given class, N + 1 states. We define a zero-state in which 
all mod .les of this given class are operable, a one-state in which one module is 
no longer working, and so on through the N-state which Indicates all modules of 
this class have failed. It is possible to calculate, from original occupation 
numbers of each of the states, what the occupation numbers for each state at any 
time in the future will be. If one starts with a zero-state completely occupied - 
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that is, probability of being in zero state is equal to one and probability for each 
of the other states up through the N-state is zero * and waits long enough, the prob- 
ability approaches zero for all the states except the N-state, which approaches one. 

The probability that, at time t, one will be able to assemble from the N components 
of a given class a set of M < N modules to perform a given task or set of tasks at 
that time is simply the probability that no units have failed, plus the probability 
that one unit has failed, plus ... and so on through the probability of a maximum 
number of failures which still provides the required hardware. It is possible then, 
by forming the sums, to determine the probability of having the necessary modules 
of a given class. If we further assume that the failures of modules in different 
classes are Independent, it is possible to write down the probability that one can 
assemble at any point in time, a complete set of hardware required by the task load 
at that time, as a product of sums. It is necessary only to propagate individual 
class probabilities, the number of which is linear with the total number of modules 

modules of type k and are 

required for the task load being considered at time t, then the w>rk can be performed 
only if < 

total available hardware. The probability that the necessary hardware is available 
is 

P(0»tt S , (1) 

k i»o, * 

where is the probability that exactly i of the type k modules have failed by 

time t. The individual P£jj(t)'s may be calculated as shown in the following derivation. 

Assume Nj^ indistinguishable modules of type k with failure rate during 
the time Interval tj^ < t < t 2 . Define as the probability that i of the modules 
have failed (and only i). Then 

for all Os i s Nj^ for given t 

( 2 ) 


and 


0 s P^^(t) s 1 


N 

L p£k(t) = 
i-0 


N, for all k. 
k 


The sets 


■^Mj^Ct)!" and define the requirements and 


in the system, rather than 2 . If there are 


1 for given t . 



Suppose that all Nj^ + 1 of the ar® known at time t^^. Any is then calculable 

at time t^ > t^ since the probability of being In the given state at t^ is simply a 
sum of products. The sum runs over all states j which the occupants of state i at 
t^ could possibly have occupied at t^ (i.e., all J < i). The products are the prob- 
abilities that the j state degenerates to the i state during the interval < t < t^, 
times the probability that the set of modules was in state j at time t^. 

Given that the set was in state j s: 1 at t^, the probability that 
it is in state i at t^ is 

cJJ'J ‘ fl - e"^k^*^2 " 


where 


The 


(N - 1) ' 

“ (Y " " jyt ' -^y-2 ■ , the biiiomial coefficient, 

probability that the set was in state j at t^^ is thus 

^ osjsisN ^ ^ 


[l - " *^1^ ] 


Finally then 


V"k<'> 


P(t) - It s 

k 1.0 




V 2 P.,(t )<k-J 

“ '■k'’- 


[t - .-Xk<^ - 


(3) 


(-'•) 


(5) 


where P(t) is the probability of having the required hardware, defined by 
}Mj^(t)|, at time t. It is assumed that the complete set of P^j^(t^)'s is 


known at some earlier time t^. 



15 



The problem of computing the reliability of performing a mission (that is, 
of successfully completing all tasks necessary for the successful completion of the 
mission) reduces to two parts. One part consists of finding the definition, for all 
points in time, of the minimal hardware set which will accomplish the task load 
required at that time. The second part consists of calculating the probability of 
having the necessary hardware available when it is required. To illustrate how one 
might use some of these concepts and techniques, it is perhaps best to describe the 
input and the output of a computer prc'rram which accomplishes the first step. 

The input to the program whicn calculates the hardware requirements as a 
function of time during the mission consists of three parts. The first part is a 
listing of what hardware is available, how many copies of each type, its power and 
weight, and initial reliability and the performance characteristics of that particular 
kind of module* In order to be as general as necessary for the treatment of many 
different kinds of systems, all modules are characterized by a set of four performance 
parameters. These are: (1) receiving rate, (2) processing rate, (3) storage 

capacity, and (4) transmitting rate. For some modules, one or more of these 
parameters may be inapplicable and, therefore, ignored by the program. A memory, 
for example, has a receiving rate, a storage capacity, and a transmittal rate, but 
no processing rate. A data bus has a receiving and transmitting rate (the same, of 
course) but no storage capacity and no processing rate. The processors have receiving, 
processing, and transmitting rates but no storage capacity. There follows, in the 
input, a list of tasks which may be required at some point in the mission. Tasks 
normally have certain hardware requirements which include, not just a list of hardware, 
but also performance characteristics required for the hardware. For example, a 
given task might have a certain memory requirement, or a certain minimum transmittal 
rate for data. These numbers are associated with the task together with the required 
modules. The program, in associating requirements with equipment actually available, 
allocates resources to best fulfill the needs of the mission. The third section of 
data names the tasks and when they are to be initiated and completed. The output 
warns of any situation which may arise in which the tasks called for cannot be 
performed due to insufficient hardware and will inform the user of excess hardware 
(if any) which is carried and the consequences in terms of excess weight and 
power requirements. The program reads che tasks sequentially and assigns hardware 
based on the requirements of each task. The algorithm for the allocation of 
hardware is somewhat arbitrary and must be designed for whatever executive system 
is being used. The algorithm employed for the sample calculation adopts the 
following guidelines: 
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(1) I£ a taak requires more than one of a given module 
type, only currently unassigned units may be con** 
sidered. If a sufficient number of such units 
can be found to accommodate the task requirements, 
the load is divided equally among such units. 

Considering only unasslgned modules minimizes the 
number of units allocated to a particular task. 

(2) If only one of a given module type is required 
by a task, then no more than one of a given type 
can be assigned to that task. A task may not be 
allocated, for example, ten percent of each of 
two memory modules. 

(3) All tasks are allocated the necessary hardware 
in the order in which the tasks are specified in 
the data deck, subject always to the restrictions 
given by (1) and (2) above. 

Sample Calculation 

Figure 2 shows the suite of ux>dules chosen for this exercise. Along 
with the module names and unit designators (A, B, C, ... etc.) the active and 

dormant failure rates, the four performance specifications previously discussed, 
and weight ana power are given. 

The hardware specifications given do not correspond to any real, specific 
hardware. Since the purpose of this example is merely to demonstrate a capability 
rather than to analyze a given hardware suite, all numbers are reported in 
arbitrary units rather than pounds, hours, watts, bits per second, etc. 
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Figure 3 llluscrates Che demands on Che hardware. For each task, 

Chere is a Use of Che necessary modules and Che performance require':'. In Che 
firsc cask, CiCled Receive Ground UpdaCe (Reev Gnd UpdC), Che CelemeCry is 
required co accepC 400 uniCs of informacion per unic Cime. The On-line Memory 
musC score 64000 unics of informacion. 

Figure 4 shows Che mission load (chac is, Che number of modules of 
each Cype required ac discrete times in the mission). Ac these times, ic also 
shows Che failure state probabilities, the Equation 5. Ac 4999. bO time 

units into Che mission Che Control System Monitor (Control Sys.-Mon) Cask is 
initiated. 

The system requires a single Inertial Measurement Unit (IMU). Since two 
such units are onboard, there are three possible states: (1) both are still working 

(probability ■ 0.6065); (2) one has failed (probability * 0.3445); or (r; both 
have failed (probability * 0.0489). The probability Chat the need for an '.MU 
can be satisfied at 4999.90 time units into Che mission is 0.9510 (0.6065 + 

0.3445). 

Figure 5 provides a summary of the utilization of available modules 
in the cas<s of no failures. It provides no direct reliability information but 
is useful iu planning. 
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FIGURE 3. COMPONENT REQUIREMENTS OF INDIVIDUAL TASKS 
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FIGURE 4. BQUIWENT ALLOCATION AND RELIABILITY STATE PROPAGATION 

(Continued) 
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T<ak 3 Invgtlaatc Scout Aatrlonlea RMulramenti 


The current Scout launch vehicle utilizea an open loop autopilot 
for guidance during lower stage actions and a spinning fourth stage. Without 
a closed loop navigation system, Scout Injection errors are considerably 
larger than for other current launch vehicles. 

Over recent years, VSPRE has Investigated a number of possible 
modifications to the Scout guidance (References 2 and 3). These modifications 
have ranged from an onboard closed loop inertial navigation system as part 
of the proposed new vehicle (ASLV) to relatively simple modifications such 
as a velocity meter and velocity vernier and/or active steering of the 
spinning fourth stage. 

The desirability of any of these modifications hinged on the 
mission requirements for accuracy. To properly assess the situation, 
accuracy requirements must be known for both currently planned launches, 
and even more Importantly, any missions which might be reassigned to 
Scout from a higher peformance vehicle. As explained In more detail under 
Task 4, adequate mission requirement data were not available for this last 
purpose. This, coupled with constraints on funding for launch vehicle 
modification, have left the decision regarding major Improvements to 
Scout accuracy open at this time. 

Spun Stage Simulation and Error Model 

In the course of Investigating alternative guidance approaches 
for the spun fourth stage, a spinning stage simulation, and a spun stage 
linear error model were developed. The simulation Involves a separate 
FORTRAN program that considers only the fourth stage. The linear error 
model Is compatible with EOMP and can easily be added to EOMP when the 
need arises (Reference 3). 

Investigation of the Dynamic Behavior of the Red Shift Probe 

VSPRE was directed by the COR to aid MSFC in assessing certain 
potentially significant error sources of the spun Scout fourth stage when 
applied to the Gravitational Red Shift Probe. 
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The Gravitational Red shift payload Is scheduled to be launched by a 
Scout vehicle In 1975. The payloa'* will remain attached to the burned»out fourth 
stage of the Scout and Is supposed to rotate at about 10 radians per second while 
describing a ballistic trajectory with a coast time out of the atmosphere of 
approximately 10,000 seconds [Reference 19]. Because of the long coast time, 
effects that are usually Ignored for short-life missions required Investigation. 

BCL was requested to Investigate the effects of several potential disturbing forces 
on dynamic stability and angular rate of the Red shift payload attached to the burned- 
out Scout fourth stage. The error sources Investigated Included 

(1) Eddy currents induced by Earth’s magnetic field 

(2) Venting of ansnonla front the side of the payload 

(3) Fourth-stage tlpoff 

(4) Fourth-stage thrust misalignment 

(5) Fourth-stage burn-out mass unbalance 

(6) Payload static and dynamic unbalance. 

The eddy current problem was handled analytically; the remaining error sources were 
Investigated by using a computer simulation of the Scout fourth stage with the Red 
shift payload attached. 

Eddy Current Damping 

Because of the long ballistic coast time, the effect of spin *ate damping 
due to eddy currents Induced by the Earth's magnetic field was Investigated. 

Reference 18 gives the following formula for the damping torque Induced In a cylinder 
spinning about Its axis of synanetry: 
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dyne*>cm 


where; 


J '* magnetic field component perpendicular to the spin axle 

a 

In Oerateda 

f ■ flneneaa ratio (length/diameter) 

(I) “ apin rate in radiana per aecond 

V ** conductivity of vehicle akin in abmi'.os/cm 

T •• akin thickneaa in cm 

r ■ radiua of the cylinder in cm. 


( 6 ) 


The following formula for magnetic field intenaity ia given in Reference 19 and 
deacribea the Earth* a magnetic field within about 10 percent accuracy: 
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where; 


J - 

o 

- 

A ■ 


horl 2 ontal intensity at the Equator ■ 0.31 Oersted 
radius of the Earth 

distance from the center of the Earth to the point 
of interest 
polar angle. 


When Equation (7) is substituted inti.' Equation (6), the result is 
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( 8 ) 


when 9 is assumed to be 90 degrees. 

The model of the Red shift payload used with Equation f8) was takex from 
drawings of the payload and the Scout fourth-stage Information from References 20 and 

21. It was assumed chat the combined package «ras an aluminum cylinder 17 inches 
in diameter and 104 inches long, spinning at 100 revolutions per minute. Skin 
thickness, t, was chosen conservatively at one-half centimeter. Using these 
numbers. Equation (8) reduces to 


T 



Ib-ft. 


(9) 
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According to Reference 22, the coast starts at 196.50 seconds at a 
geocentric radius o£ 22.75 X 10^ feet (with Earth's radius being 20.9 X 10^ 
feet). At this altitude, ■ 0.61, and T • 6. A X 10 ^ Ib-ft. 

Decay of the body's angular rate is given by 
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(I) e 
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( 10 ) 


where H is the angular momentum of the body about the spin axis. At the start 
of coast, H ■ 27.5 lb*ft-sec for the modeled Red shift payload. Figure 6 shows 
the geocentric radius (A) of the payload as a function of time. Tt also shows 
the stair-step approximation to A used incrementally in Equation (9) to obtain 
an approximation to the roll rate damping torque during coast. Equation (10) may 
be rewritten for piecewise solution as: 
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where 


and 


H 


1.05 X 10 
2.63 u) 

n 


-3 


Results of applying Equation (11) to the stair-step approximation of A are given 
in Table 3. The decay of the roll rate is plotted in Figure 7. It is readily 
seen that the decay due to eddy current damping is predicted to be less than one 
percent for the 10,000 seconds of the experiment. 

Spinning Stoge Simulation 


Previous work for VSPRE included derivation of equations of motion and 
development of a 6 degree-of-freedom simulation for spinning stages (Reference 3). 
The pertinent portions of this effort are repeated here for convenience. 










100.25 



FIGURE 7. SPIN DECAY DUE TO EDDY CURRENTS 
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Equations of Motion 

The general equations of rotational motion^ < presented in 
Reference 23, are expressed for our notation in Equation 12: 


T « J (i) + S(uj) Jo} + 


where 

T = the applied torque vector, 

J = the moment of inertia about the center of 

gravity, a 3 by 3 matrix, 

0 ) =5 the vector of angular rates, 

S( ) = the skew symmetric operator, a 3 by 3 
>v 

matrix, 

= the ith engine flow rate, and 

r . ~ the location vector of the ith engine, 

ei 

relative to the c.g. 

Equation 12 can be rewritten as a differential equation for the angular 


rates: 


0) = J ^ ^ T - S(oj) J 0) - D uo J 


The thrte-by-three matrix D is used to represent the last two tenns 
of Equation 12: 


j+E*i 


^ S(y)x = y cross x, the vector cross product. 


/ 
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Equation 14 has terms for the time rate of change of the moment of inertia, 

• 

J, and the angular momentum loss due to mass expulsion, jet damping. As 
written, Equation 14 assumes nozzles of zero area. Figure 8a represents a 
stage with four nozzles of zero area, each located away from the roll axis. 

If the nozzles are not on the roll axis, the jet damping effect is significant, 
even for zero area nozzles. For a more practical development, especially for 
the case of one nozzle on the roll axis, represented by Figure 8b, an integral 
over the exit area of the nozzle(s) is required. Consider a nozrle of area A, 
located with its centroid at 1 relative to the vehicle c.g. as shown In Figure 9. 
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FIGURE 8a. SEVERAL NOZZLES OF 
INFINITESIMAL AREA 


FIGURE 8b. ONE LARGE NOZZLE 
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FIGURE 9. INTEGRATION OVER NOZZLE AREA 
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Let p be the mass flow per unit area, and s the vector from the nozzle 
centroid to an element of area, dA: 

*1 " J* P + 8) dA (15) 

A 

Expanding the skew squared gives 

*1 ®^<^el^ • I P [ S(s)S(f) + S^s) ] dA ( 16 ) 

If the mass flow density, p, Is Independent of the location, s, the 
first term of Equation 11 can be expressed in terms of the total mass 
flow rate ib. Since s Is relative to the centroid of the area, 

f S(s)dA = 0 (17) 

A 

and the two middle terms are zero. The last term is a function of 
the planar moment of the area. Therefore, 

D = i + S^(sp) ( 18 ) 

1 

2 

where is the location of the ith nozzle and S (s^) is a function of 
nozzle area and orientation. 

Six Degree of Freedom Simulation 


To perform analyses , and to provide a basis for 
evaluating linear analysis assumptions, a general spun stage 
six-degree-of- freedom (6 DOF) simulation was written, xhe simulation 
performs a fourth-order Runge-Kutta-Gill integration of the 
position, velocity and angular rate equations: 

t = V, 

= B f/m, and (19) 

(b = J -S(u)) Ju) - j(U ^ ^ 1 

1 

where r and v are inertial position and velocity vectors, B is the 
three by three direction cosine transformation from body to 
inertial space, f is the applied force in body coordinates, and m is 
the vehicle mass. The attitude equation, 

B » BS(ui), 
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Is not integrated. A closed form solution, developed for use In error anelytts» 
Is used. This solution Is valid over any time Interval for periods of 
simultaneous constant angular rates. While only an approximation for 
this problem (due to time-varying rates), the error is extremely small. 

Its use Insures maintenance of the orthogonality of B, and its execution 
efficiency is excellent. 

The time-varying coefficient arrays are computed from data 
describing the vehicle on a piece part basis. Each piece part Is 
considered a hollow cylinder with its axis parallel to the roll axis. 

For each of these cylinders, the following data must be provided: 

1 . Length 

2. Inside diameter 

3. Outside diameter 

4. Initial weight 

5,6,7. Location of cylinder c.g. in body coordinates 

8. Ignition time 

9. Mass flow rate 

10. Nozzle exit area, and 

11,12,13. Thrust vector in body coordinates. 


For inert masses, items 8 through 13 are not used. Fuel 
masses are assumed to burn from the inside out. This approximates the 
burning of many solid fuel engines. The simulation computes the 
system moment of inertia, c.g. location, damping terms, etc., by summing 
whe appropriate data for each piece part. Analysis of a Scout fourth stage 
and its payload, using this cylindrical model, produces parameters that 
agree well with published values (Reference 24), 

Information used in simulating the Red shift payload was taken from 
References 25 - 28. Mass property data were contained in Reference 25, but these 
data were modified. The weight given in this reference was 144.1 lb; late 
information gave the weight as 165 lb. The moments of inertia given in 
Reference 25 were increased by the ratio • 


Any assumed burning pattern can be simulated with minor coding changes. 
** Telephone conversation with Mr. Ernest Nathan of MSFC. 
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The Scout fourth stage was modeled after information taken from References 20, 

21 and 29. Ihe model of the stage consisted of four pieces: 

1. Engine case 

2. Fuel 

3. Nozzle 

4. Upper D section. 

The parameters used to describe these pieces and the payload are shown In Table 4. 


TABLE 4. VALUES USED IN SIMULATION 



Payload 

Engine 

Nozzle 

Fuel 

Upper D Section 

Length (In.) 

35.0 

48.0 

34.0 

36.7 

35.2 

Inside Dla. (In.) 

4.0 

8.4 

8.4 

3.0 

9.6 

Outside Dla. (In.) 

8.0 

8.6 

8.6 

9.04 

10.8 

Initial Wt. (lb.) 

165.0 

30.0 

21.2 

613.9 

16.48 

eg X (in.) 

26.9 

58.3 

89.0 

67.18 

91.33 

eg Y (In.) 

0.0 

0.0 

0.0 

0.0 

0.0 

eg 2 (in.) 

0.0 

0.0 

0.0 

0.0 

0.0 

Ignition time (sec.) 




0.0 


Mass flow rate J 

\sec / 




-21.2 


Nozzle exit area (ft^) 




2.25 


Thrust X (lb.) 




5925.0 


Y (lb.) 




0.0 


Z (lb.) 




0.0 
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The simulation program Is written to provide graphic output. An 
example Is shown In Figure 10. The Y-Z plane Is defined as being normal to 
the vehicle roll axis at the start of the simulation. It then translates with 
the vehicle center~of~gravlty. The trace on the Y-Z plane, then, is a function 
of the angular motion of the vehicle's roll axis, as shown In Figure 11. The 
dimensions on the Y and Z axes correspond to the sines of the angles that the 
roll axis makes with respect to the Y and Z directions on the Y-Z plane. 

Discrete events are shown as triangles on the output. Each plot has 
a triangle at the origin corresponding to Ignition of the fourth-stage engine. 

In order to avoid putting too much data on a single plot, a number of plots may 
be used, each covering a distinct time period. The triangle at the origin Is 
shown for reference on all plots. Including those that start at a time greater 
than that of fourth-stage Ignition. When a multiple-segment burn Is used to 
approximate the engine, a triangle Is printed at the start of each segment's 
burn. 

Although results obtained by the simulation appear to be reasonable 
during the powered flight portion of the mission, an unwelcome anomaly was found 
to occur at burnout. The pattern in Figure 12 is typical for Scout fourth stage with 
thrust misalignment during powered flight; however, the "tail" seen In the lower 
left of the figure Is definitely not to be expected. The pattern should be a 
circle with diameter about .068 (or about 4 degrees) after termination of powered 
flight. This discontinuity Is apparently due to assuming a rectangular burn 
with constant thrust and flow rate for the burn interval. When an eight-segment 
staircase approximation to the burn was used, this effect was minimized. Further 
discussion of this modification and Its results is contained later In this 
section. 

In order to obtain appropriate tipoff response, it was necessary to 
modify the fourth- stage exit area to 2.25 square feet rather than the 1.5 square 
feet given in Reference 20. 
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Error Sources Conildered 

Several error eourcea were ilnulated in order to determine their 
effect on dynamic behavior of the Red shift payload. These sources were 

(1) Ammonia venting from the payload 

(2) Scout fourth-stage tipoff 

(3) Fourth-stage thrust misalignment 

(4) Fourth-stage burnout unbalance 

(5) Payload static and dynamic unbalances. 


ii. 


f i 


Ammonia Venting. The payload includes a supply of liquid wnmonia for use 
in cooling the maser cavity. After the ammonia changes to gas during the cooling 
process it is vented overboard, normal o the roll axis, at Scout Station No. 29.97. 
Gaseous ammonia is vented at a rate of 0.28 pound-per-hour with a resultant force of 
0.0034 pound (Reference 30). Figure 13 shows the results of the ammonia venting 
for the first 100 seconds of operation following fourth-stage ignition. It is 
seen that the original disturbance which results in a half-cone angle of about 
0.6 X 10 degree, decays by the end of fourth-stage thrusting to about 0.2 x 10*'’ 
degree. During the coast, this angle is maintained, although modulated by the normal 
force due to the ammonia venting. 


u 

; 5 


Fourth-Stage Tipoff. Figure 14 shows the eTfect of fourth-stage tipoff. 

The initial response is a half-cone angle of about S degrees, which is damped to 
about 1.1 degrees by the time of burnout. This corresponds quite well with the values 
of 5.22 degrees at ignition, given in Reference 31, and 1.11 degrees at burnout, given 
in Reference 29. 


Thrust Misalignment. The effective 3-cr thrust misalignment for the FW-4 
engine is given as 0.3 degree in Reference 29. When this value is included in the 
simulation, the plot of Figure IS results. The initial half-cone angle of about 4.6 
degrees decays to about 1.8 degrees at engine burnout. The "stem" going to the lower 
left-hand comer of the plot Is a result of simulating the engine bum as a single, 
rectangular bum. A more detailed simulation would show the 1.8 degree half-cone 
angle continuing during coast. 


Z-KXIQ -3) 



I 

i J 

^ \ 


OI 

o 


.019 


-T— 

-.014 


“1 

-.009 


I » \ 

-.005 -.000 .006 

Y-UXIO -3) 


.009 


.014 


■019 


FIGURE 13. EFFECT OF AMMONIA VENTING FOR FIRST 100 

SECONDS AFTER FOURTH-STAGE IGNITION ^ 
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Fourth»Stage Burnout Lnbalance. This error source Is attributed to 
blistering of the case insulation (Reference 29), Although a stnall amount of flight 
data has resulted in a 3-? dynamic unbalance value of 2400 oz-in^ for this source, 
the LTV design criterion calls for an upper limit of 8000 oz-in2. This latter value 
was used in the simulation. The effect of this error source was found to be inslg~ 
nif leant and no plot is Included in this report. 

Payload Unbalances. Reference 20 gives 12 oz-in and 200 oz-in for 
maximum values of static and dynamic unbalances, respectively. When these terms were 
added to the simulation, no significant error was found to result. Therefore, no plot 
is included of the effect of these error sources. 
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Thrust Misalignment and Tipoff. Since fourth-stage thrust misalignment . ! 

and tipoff were the major error contributions, runs were made with both of these error 
sources combined in such a way as to produce maximum error. During these runs it was 
found that the initial half-cone angle was about 8.5 degrees as can be seen from 
Figure 16. 1 - 

During this series of runs, it was decided toinvestigate the cause of 
the "tail" following engine cutoff. Three-segment and eight-segment models of the . ; 

thrust and flow rate for the engine were Incorporated into the simulation (See Figure 17-'^> 
Results for one-, three-, and eight-segment models during the last portion of powered 
flight and the first few seconds afterward are shown in Figure 18-20, respectively. For t 
the original, one-segment model (Figure 18), a half-cone angle of 20 degrees results. .1 
Figure 19 shows a half-cone angle of about 9.5 degrees for the three-segment model. 

I 

When the 8-segment model is used (Figure 20), the half-cone angle is about 2 degrees* 

This is the same value that wuld be predicted by examining the trace of Figure 16 just ; 
before burnout* Thus, it is concluded that the anomalous behavior in many of the plots 
after burnout can be ignored* 
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Concluding Remarks 

From the results of the Investigations, it appears that the error 
sources considered are not likely to cause significant problems in conducting 
the Redshlft mission. It must be realized, however, that the simulation used 
for most of these studies does not use detailed models of the vehicle and 
payload and, thus, subtle effects could remain undiscovered. 
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Task 4 Perform Mission Requirements AaBessments 

Studies conducted under this contract have been concerned with 
the assessment of conventional launch vehicle astrlonlcs capabilities and 
their adequacy for projected mission requirements. Injection accuracy has 
been a primary concern. Navigation accuracy analysis techniques for both open 
loop (Scout, Titan HID), and closed loop (Delta, Atlas/Centaur, Titan/ 

Centaur, etc.) were developed. (References 1-3) 

These techniques permit evaluation of current (or proposed 
Improvements to) Injection accuracy. To objectively apply these techniques, 
a concise set of allowable injection conditions must be known. Investigations 
of missions frequently indicate arbitrary or vague specification of allowable 
deviation from the desired aim point. For example, one payload was placed In 
an orbit well outside the specified "window" but was found to provide experi- 
menters with more data than would have been obtained had the nominal orbit 
been achieved. Considerations such as these led to the Investigation of 
the relationship of experiment oriented parameters such as orblt/earth and 
orbit/sun precession rates to satellite injection conditions. 

Concurrently, Battelle personnel were Involved to varying degrees In 
the Initial planning phases of future missions (HEAD, SEASAT and LAGEOS). In 
each of these cases, potential experimenters needed answers to questions 
requiring analyses outside their Individual technical disciplines. Thus, NASA 
launch vehicle offices and vehicle prime contractors were called upon to provide 
analytical support well In advanced of the approval of experiments and missions. 
Conducting effective tradeoff analyses of even a few alternative orbits can r 
require significant manpower and computer time If traditional detailed analytical 
techniques are used. 

To meet both the needs of the astrlonlc evaluation function and the 
mission planning function, the Interactive Graphics Orbit Selection (IGOS) 
program was written. 

Experience has shown the value of presenting a synopsis of system 
propulsion capabilities and limitations to potential space transportation 
system users. This was accomplished by periodic publication of. Reference 34, 
the Launch Vehicle Estimating Factor Book (EFB) by the OSS Launch Vehicles 
and Propulsion Programs Division, Code SV [now Expendable Vehicle Programs 
(Code MV) in the Office of Space Flight (OSF)] • This planning handbook (EFB) 
contains data, usually In graphical form, for: 
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(1) Launch Vehicle payload capabilities versus 
characteristic velocity, 

(2) Earth orbit and planetary mission Vq requirements 

(3) Determining Launch Vehicle payload capability versus 
Earth orbit parameters (perigee, apogee, and Inclina* 
tion) for specific vehicles 

(4) Mission restrictions (orbit life time, range safety, 
etc.). 

These data are presented in a form suitable for Interpretation 
by a payload planner. Vfhlle this necess'l tates reasonable engineering 
approximations, the data are sufficiently accurate for trades between the 
users desires and the transportation system capabilities. Developing mission 
specifications using this handbook provided an opportunity for the payload 
planner to consider a wider range of alternatives than would be possible using 
detailed analysis techniques. Detailed analysis could be performed later In 
the planning cycle and be based on a feasible set of mission specifications. 

The availability of Interactive graphics terminals and central com- 
puter software has permitted developing various planning tools Including one 
for preliminary mission planning which Is far more effective than the tradi- 
tional handbook format. Several years of successful operation of the NASA 
Interactive Planning System (NIPS) by both OSS and OA had previously demon- 
stratedthe feasibility and reliability of financial analyses performed us^ng 
Interactive graphics terminals for decision making, report generation, and 
data file maintenance. The interactive graphics tools have now been extended 
to permit preliminary mission planning by this program Is far more versatile 
and effective for Earth orbital mission planning than the traditional handbook 
format. 

A detailed users manual and sample work session are Included in 
Appendices A and B. The following paragraphs summarize the operation of 


IGOS 
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IGOS Program Description 

The IGOS program faclllcatea qulck-response asaessment of Earth> 
orbit mission requirements and their trade<*o££s with launch-vehicle capabilities. 
By using an interactive graphics computer terminal, a quick-response low-cost 
analysis can be initiated early in the mission planning cycle. Application 
of this program in no way eliminates the need for the traditional detailed 
mission analyses during later phases of mission planning. It does, however, 
reduce the need f.r>r major iterations in the planning cycle due to unacceptable 
orbit characteristics. 

IGOS utilizes the capabilities of interactive graphics terminals, so 
experimenters and mission planners can clearly visualize the many possible 
orbit alternatives that can satisfy their mission requirements. The operating 
environment for the IGOS system is Illustrated in Figure 21. The user 
enters a set of mission requirements at a remote graphics terminal and 
receives appropriate graphical and textual displays. The information 
displayed enables the user to select acceptable alternative Earth orbits. 

Graphical output is in a two-dimensional design space - altitude and 
inclination for the vertical and horizontal axes. Thus, a point in the plot 
represents a circular orbit. Regions representing unacceptable orbits are 
shaded, as shown in Figure 22. The shaded areas result from the constraints 
imposed by the following requirements and phenomena: 

(1) Earth -observation coverage 

(2) Radiation environment 

(3) Sun-orbit precession 

(4) Orbit decay 




FIGURE 21. lOOS OPERATION 












FIGURE 22, TYPICAL IGOS GRAPHICAL OUTPUT 
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(5) Launch**vehicle injection accuracy 

(6) Launch'vehic le payload capabilities. 

Each of these is described in the following paragraphs. The 
user-supplied data are listed and a sample plot is shown. This is followed 
by a hypothetical exercise using all of the features. 

Earth-Observation Coverage 

Input Data: 

e Latitude range of intet'est, 
e Days allowed for ^^iewing all longitudes 
within the latitude range, and 
e Sensor data (maximum slant range, minimum 
elevation above the horizon, field of view, 
and resolution). 

In the resulting display, the shaded areas indicate all orbits 
which do not view every longitude in the latitude range of interest in the 
allowed time. 

Radiation Environment 


Input Data: 

• Aluminum shielding density, 
e Allowable fluence (flux-time integral) 
accumulated inside the shielding, and 
e Mission duration for .accumulating the fluence. 
Orbits which experience excessive fluence are indicated by 


verti':al hatching in Figure 22. 
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Orbit to Sun Precession 

Input Data: 

• Allowable precession rate. 

The acceptable region is shown by dashed vertical lines, as shown 
in Figure 23. The center line denotes the locus of zero precession (sun 
synchronous orbits), and the other two lines indicate the allowable limits 
for an acceptable precession rate. 

Orbit Decay and Launch-Vehicle Injection Accuracy 

These separate phenomena are both related to a specific nominal 
injection condition, and their effects are displayed together. 

Input Data: 

• Injection altitude, 

• Injection inclination, 

• Spacecraft ballistic coefficient, 

• Launch date, and 

• Mission duration. 

A vertical line is then displayed at the specified inclination, as 
shown in Figure 24. "Tick** marks indicate the nominal, plus and minus, one- 
sigma and three-sigma altitudes. Ticks on the left represent the dispersion 
of initial orbits; ticks on the right represent the corresponding altitudes 
after the specified decay time. Ticks occurring above the top of the plot 
are indicated by displaying their altitude; those below the bottom are 
annotated with the date they would pass the minimum altitude. 
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Launch-Vehicle Payload Capabilities 

Input Data: 

• Launch vehicle and 

• Payload weights of interest. 

All launch-vehicle and launch-site data are stored with the program. 
When a launch vehicle is selected, equipayload weight contours are plotted, as 
shown in Figure 25. The contours are the altitude at which the vehicle can inject 
the payload into a circular orbit, as a function of inclination. For each 
inclination, all of the possible launch sites for the selected vehicle are 
considered, and the site yielding the highest altitude is chosen. Launch azimuth 
c' straints and plane change losses are included in the performance calculations. 

Although Figures 22 through 25 show separately the important features 
of IGOS, the full benefit of the program is apparent when these plots are 
superimposed and modified interactively. Consider the following problem. 

It is desired to view all the Earth's surface between 20° and 60°N 
latitude at least once every 30 days for a 200-day mission. Assume a sensor 
with the following characteristics: 

• Maximum slant range 1200 n.mi. 

• Minimum elevation 40° 

• Field of view 40° 

• Resolution of 10 n.mi. on the Earth * 0.1° at the 
sensor. 

The orbit should precess relative to the sun less than 30° during 
the 200-day mission. 



> 4 ^ 
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2 

It is feasible to carry 0.15 g/cm of shielding to protect critical 

12 

subsystems, and it is desired to keep the fluence below 1.24 x 10 1-Mev 

2 

equivalent electrons per cm over the 200-day period. 

A first look at these requirements results in the plot shown 

Figure 26. Study of the display indicates that orbits of interest are near 
90^ inclination. Figure 27 was created to **zoom in** on those inclinations. 

A trial nominal orbit of 185 n.mi. altitude, 96® inclination, and a 
tentative launch vehicle (hypothetical data are shown) are selected, which 
result in the solid vertical line showing injection dispersions and orbit decay 
being added to the display (Figure 27). That line shows that the expected 
injection dispersions lie across two unacceptable regions (insufficient coverage). 
Also, there is noticeable decay in 200-days, especially for a -3a injection. 

Using the indicated decay as an estimate of decay rate, an injection near +2a 
would pass through the narrow unacceptable region. Likewise, injection near 
-3a would pass through the larger unacceptable region. While this region is 
wider, the decay rate at the lower altitude is faster. More detailed evaluation 

of this problem would be possible by plotting the decay for shorter mission 
times to establish the time spent in the unacceptable regions. Likewise, 
other coverage requirements could be analyzed to establish what coverage is 
achieved in the unacceptable regions. 

This example demonstrates the capability of IGOS to serve as an 
automated sketch pad for rapid assessment of alternative orbits. The avail- 
ability of such a tool offers several benefits. First, the quick-response 
analysis will ensure that a wider range of alternatives is considered early 
in the planning cycle, thus the orbit finally selected will be more likely 
to be a best match to the all of the objectives of the experimenters. Second, 
the improved preliminary planning cycle will reduce costly Iterations during 
the later detailed analysis phases of mission planning. 




t 

V 


* ■ 






FIGURE 26. SAMPLE PROBLEM COVERAGE, RADIATION, AND SUN-ORBIT PRECESSION REQUIREMENTS 






FIGURE 27. SAMPLE PROBLEM WITH TRIAL ORBIT AND LAUNCH VEHICLE 
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This program displays a concise set of mission requirements 
in terms of allowable injection dispersions. This format provides the 
mission specifications necessary to assess the suitability of booster 
guidance accuracy for both current systems and proposed modifications. The 
program is currently operational on the Battelle-Columbus computers with 
terminals in Columbus and at NASA Headquarters and is on computers at the 
Marshall Space Flight Center. 

IGOS Mathematical Models 

IGOS uses a set of separate mathematical models, each with its 
o\>m input parnTnct''rs .ind uach computing data for plotting on a common 
set oC scales. The urogram is organi^ed into groups of subroutines, each 
associated with a separate mission analysis phenomena. These groups of 
routines can be Lnought of as separate analysis programs. Each generates 
data for plotting in the following two-dimensional design space: 

Altitude h, h . h < h , and 
mxp ~ - max 

Inc lint ion, [, T . f I 

min •“ -nax 

The following discussion describes each of the separate models. 


Ea r th (.overage Model 


Hie object of the Earth coverage model is to shade those portions 
of the (Icbign space (I and h) which do not satisfy the requirement to view 
all points on the Earth’s surface between latitudes 6^ and 0^ least once 
in any time period, T^. 

The analysis considers a sensor model which computes the radius 
of the circular area on the Earth’s surface viewed at any instant in time* 
fliis radius is expressed in ti rms of the angle, ' , of a cone from the 
Earth’s center and is a function >f chc orbit altitude and the sensor 
ci>a rac tc r i s t ics as discussed in the Sensor Model section later in this 
report . 

Wlien a single latitude is ansidcred, the swath viewed by a 
s 'nsor iniersects the latitude ir ceu or U o hands of Iracitude with 
width V. V is a function of the englt , the latitude and the orbit 
iiltitude and incl i e.ation. 








¥ 




k 


70 

The ground crack o£ Che orbic shifCs durlag each orbical pass 
due to Che EarCh's rocacional race, lu , and Che orbic 's precession race, 
u)^. Each pass is shifted in longitude by an angle, S. After many orbital 
passes, a distribution of c^^ound track/latitude intersections resulcs. 

The spacing between these intersections are referred Co as gaps, G. If 
the largest gap is less chan the longitude viewed, total coverage is 
achieved. T'.aC is, when 

G (T , h, I) < V ti|((h), 0, h, I] (20) 

M c 

is satisfied, adequate coverage is achieved. 

The following sections discuss the sensor model with viewing 
angle the longitude viewed (V), the orbit shift (S), and the maximum 

gap (G ) in that order. 

M 


j 




i. 


Earth-Observation Sensor Model. IGOS includes a sensor model 
which is used to calculate the instantaneous area viewable as a function 
of satellite altitude. Among the parameters included are the following: 


(1) 

0 * 
0 

sensor field of view limit 


(2) 

a = 
0 

sensor angular resolution 


(3) 

s = 
o 

maximum sensor-to-target slaut range 
(may be limited by uplink sensitivity or 
downlink power supply) 


(4) 

T i, 1 1 n 

Y =* 

minimum elevation angle (must be set to 
meet minimum lighting constraints, to 
avoid excess refraction of light signal, 
and to clear ground obstructions). 



Items 1 and 2 are sensor parameters. Items 3 and 4 may be either sensor 
parameters or mission design factors* 


/ 






k 




V 
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Other 

following: 

variables and parameters used in the model Include 

(5) 

s 

■ slant range 

(6) 

h 

* altitude 

(7) 


* viewed area half -angle 

(8) 

Y 

“ elevation 

(9) 

e 

* TT - 2 (t + Y) 

(10) 


® linear resolution at target center 

(11) 

X. 

* linear resolution at target perimeter 

(12) 


® earth radius 

(13) 


Lq = flags which indicate which parameter 
limits the radius of the viewed area, 


Figure 28 depicts some of the relationships among these variables. 

The constraints to be satisfied are: 

(14) s < 

(15) e < 0^ 

(16) Y > 

(17) < s A^/sin Y 

As defined here, a smaller value for resolution indicates "better" resolution. 
Generally, each restricted parameter defines t as a function of altitude. The 
composite graph of ^ji as a function of h is the piecewise-smooth curve repre- 
senting, at each altitude, the most restrictive parameters. Figure 29 shows 
the flow of the calculation and illustrates the equations used. The limiter 
flags are initialized to zero. If the altitude is greater than the maximum 
sensor-to-target slant range, the viewed-area half-angle is zero. Otherwise, 
it is assumed that the limiting parameter is the elevation restriction and 
theta is calculated. If the sensor field of view limit is exceeded, theta 
is set equal to it and the elevation angle is recalculated. In either case, 
if the resulting slant range is too great both theta and the elevation angle 
are recalculated after setting the slant-range-to-target-periraeter equal to 
its limiting value. 



To illustrate the utility of the model, an example is presented. 

The sensor is assumed to have a field of view ‘ 40 degrees and an angular 
resolution of 0.01 degree. The desired lineal resolution is 1 nautical mile. 
The minimum elevation angle is 40 degrees and the range limit is 1200 nautical 
miles. 

In Figure 30, the solid curves show i|f as a function of h considering 
each limit as If it were the only constraint. The cross-hatched segmented line 
is the piecewise-smooth curve reprasenting as a function of h. In this case, 
the resolution constraint is never a limiting factor for altitudes ranging 
up to 1000 nautical miles. The sensor fleld-of-view is limiting for the lower 
altitudes. At 660 natical miles the elevation angle is the most restrictive. 

At 870 nautical miles, the range restriction is dominant. The maximum value 
of (approximately 12 degrees at 870 nautical miles altitude) corresponds to 
a great circle distance of about 720 nautical miles on the viewed area. 
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Longitude Viewed, A single orbit pass covers the surface 
of the Earth with a swath as shown in Figure 31 « The band of longitudes 
vi^'wed, V, is shown for a particular latitude, 9 * Note that for a 
larger inclination the longitude viewed would be two smaller bands, while 
for a smaller Inclination no longitude would be viewed. 

To compute the portion of any latitude viewed in a single 
satellite pass, it will be convenient to develop the orthogonal trans- 
formation from a set of orbit related coordinates, to Earth coordinates, 

r . This transformation is 
e 
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I 

I 


( 21 ) 


where 


c. and 


the cosine and sine of the precession of the 
orbit relative to the earth, i.e.: 


Cj and Gj 
the coordinates are: 


4, = (a.p - 

= the cosine and sine of the inclination, and 


r 

-ne 



in the equatorial piano 

?n the equatorial plane, normal to x^, 

the polar axis 


and 


z 

o 


the lino of nodes, 

In the orbital plane normal to x^, 
normal to the orbital plane 
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Since all positions of interest are on the Earth^s surface, which is 
assumed to be spherical, all position vectors are treated as unit 
vectors. 

The upper and lower edges of the swath are expressed in 
orbital coordinates as 


C, c 

t o 


r « 
-0 


c, s 
^ o 




( 22 ) 


where \jf is the sensor angle, and and are the sine and cosine of 
the orbit anomaly angle, 

The edges, in Earth coordinates, are found by applying the 
transformation of Equation 21 to Equation 22: 


— — 


~ — 

X 

e 


c c c - sees ± S.StS, 

^ilfo ^Iijfo <SI<lt 

^e 


s^e.e + c.e c ,s ± c.b^s, 
(7'to ciliro (Alt 

1 


s_e,s ± e,s, 
I t o It 


(23) 


The crossing of the latitude of interest is found by setting z equivalent 

e 


to the sine of the latitude, then solving for the orbit angle 

1 

+ c^sj/s^c^ 


a) t = sin 
o 


■ [' 




(24) 


In the interval o < a^t v 2rr, Equation 24 has four, two, or no 
solutions depending on the specific angles: 

jlj < (0 - t) No solution 


(9 - *) < 

(9 + tir) V 


(9 + iO Two solutions, and 
Four solutions. 


f . 

> .i 


s 

i ^ 


, \ t 


These cases correspond to no, one, or two bands of longitudes viewed 


at a specific latitude. 
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With sets of solutions for uj^t at the crossing of the latitude, 
the corresponding longitudes can be found from 

V « tan“^ | j (^5) 

with and obtained from Equation 23. The band of longitudes viewed 
can be found by taking the appropriate differences of Equation 25 for the 
special cases of the solution to Equation 24. The resulting function is 
sketched for a single altitude in Figure 32. For at-? in IGOS, it is 
interesting to consider the value of V as a function of both altitude 
and inclination (the design space) as shown in Figure 33. 


Orbit Shift Parameter . To compute the distribution of gaps, it is 
necessary to know the amount of orbit shift per orbit period, shown as S in 
Figures 31 and 34. The shift is 


S = u) t 
pe o 


(26) 


where is the precession relative to the rotating Earth, 


U> = U) - U) 

pe p e 


(27) 


and t^ is the orbit period, 


t = 2n/(jD 
0 0 


(28) 


Combining Equations 26, 27, and 28 gives 


S = 


2n (1) - u' 


(29) 


'1 i. , ^ ^ 




V 







FIGURE 34. SUCCESSIVE GROUND TRACKS 
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The Earth's rate is a constant, a> = 7.292 x lO”^ rad/sec, while w 

e ’ p 

and u) are functions of the orbital altitude and inclination: 

0 

,1/2, 3/2 
CD * p /r 
o 


and from Reference 32 


jn 


U) = K (r /r) cos I 
p p e 


(30) 


(31) 


where 


r = orbital radius (h + r^), 

r ** Earth's equatorial radius (20902313 ft), 
e 

p. = Earth's gravitational constant (1.407528 x 10^^ ft ^^s^), and 

K = “9.97 deg/day or -2.0149 x 10 ^ rad/s. 

P 

Substituting Equations 30 and 31 into Equation 29 gives 


14 


t 


S = 2tt 




7/2 


CO 


1/2 2 

L M' r 


cos I 


e 3/2 

1/2 ^ 


M' 


(32) 


It should be noted that, for constant altitude, S is linear on the cosine 
of inclination. 


1 


S = c (r) + c, (r) cos I 
a D 


(33) 


where 


, CD 3/2 

c = 2 tt er , and 

® 1/2 


,1 

4 

I 


c - , ^ 

Cb - 2 tt P e 

1/2 2 

M. r 


I 
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Thus, for a known altitude, it Is a simple matter to cor.7«rt between 
Inclination and orslt shift. The relationship between shift and 
Inclination Is shown In Figure 35. ihe curvature of the locil of 
constant s result from the dependence on the cosine of Inclination. 

Distribution of Gaps . Figure 36 shows the subdivision of longitude 
by many orbital paases. Note that, for the first few passes, the gaps are: 

(1) 1 gap of (2rr- NS), and 

(2) N gaps of S. 

However, after a specific number of passes, the original gaps are subdivided 
Into even smaller Intervals, In general, after any number of passes, three 
gap sizes exist such that: 

^ 1^1 S ^2 

N^ + N 2 + N^ - N + 1. 

Now consider a fixed number of passes, N^. For the IGOS coverage model 

T 0) 

N 

p 2 tt 

a constant for each altitude of Interest. The coverage function requires 

calculation of the maximum gap width as a function of S [i.e., G„(S)1. 

If many cases are evaluated, each with the same N and differing 

P 

value of S, a plot such as Figure 37 will result. Note that, while there 
are generally three gap sizes, at specific values of S the smaller gap 
width goes to zero and the two larger sizes become equal. At any of 


0 


f i 






these discrete values, S^, there is only one gap size, H.ese cases 


f 








FIGURE 37, GAP SIZES VERSUS SHIFT 


are the values of S which are an exact proper fraction of 2'^, 


S 

I 



( 34 ) 


where the denominator must be less than or equal to the number of 
passes. 


The gap width is then 



(35) 


that is, with all gaps equal, the product ot width and number must 
be 2rr. These values of S are often called resonant orbits. 

With the discrete values known, the maximum gap varies 
linearly between successive points, and G^, as shown by the 
upper line of Figure 37. The fact that the maxiinwn and all other gaps 
vary linearly with S can be shown by considering any pair of longitudes. 
After 1 passes, the longitude !<* 

21 )^ = iS - 2n 


where is the largest integer < iS/2TT. is needed to adjust0^ to be in 
the interval 0,2rr. Thus, any longitude is piecewise linear on S and the 
gaps, which are merely differences (^^ - are also piecewise linear in 



To compute the discrete values, we must find all values of S 

which are exact proper fractions 'f 2 tt with denominators less than or 

equal to An example of a set of such points is shown in Figure 38. 

The points which represent proper fractions are shown solid. The improper 

fractions are shown as circles. The hyperbolae for constant M have been 

added to aio .nterpre*.?» .ion. ^ 

Fo a case with 10 passes, the maximum gap function is plotted 

in Figure 39. With 11 passes, the nlot shown with a dashed line results. 

Note thiic one additional pass greatly changes the shape of the function. 

When typical coverage times are considered, N^, can be quite 

large. For example, with near Earth orbits having periods on the order 

of 120 min (12 passes per day) and coverage times, T^, of up to 30 days, 

Np of up to 360 can be obtained. Plots such as Figure 38 and Figure 39 

would have on the order of N^I points, a very large number. However, the 

function G (S) is needed only to establish when the gap width exceeds 
M 

a known value G . To compute the S for which G (S) > G , only the point 
o mo 

pairs with at least one point greater than G^ are needed. Thus, an 

algorithm for finding G (S) can be developed which does not require 

m 

inspection of all possible points. The points will exceed G^^ when 

N. < N = G /2tt . (36) 

1 0 o 

Thus, it is necessary to find each point (integer fraction), 

S.» , with a denominator less than G^/2n. Then, the nearest adjacent 

fractions (S. and S.) with numerators less than N must be found. Using 
i i' p 

these points G (S) will result as shown in Figures 40 and 41. Wliile the 
m 

function is unknown for regions between points with less than it 
is accurate for all segments above or crossing G^. Thus, it is adequate 
for locating the values of S at any value G equal to or greater than G^. 
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A further reduction in the number of points to be calculated 

is possible because, for any one altitude, a limi':ed range of inclinations 

and, thus, a limited range of S is of interest. 

With the function G (S) known, and S determined from the altitude 
n 

and inclination, G as a function of altitude and inclination are known, 
n 

A sketch of this function is shown in Figure 42. The sketch has been 
simplified to ease interpretation. In general, the function is a series 
of ridges and valley:; which are distributed along the locus of equal orbit 
shift as shown earlier in Figure 35. 


Combined Coverage Model . The IGOS coverage model operates by 
combining the longitude viewed model (Figure 33) and the maximum gap distri* 
bution model (Figure 42) as shown in Figure 43. Adequate coverage is assured 
when V is greater than G giving a display such as that shown in Figure 44. 

The plot is generated by solving for the gap distribution (Figure 
37) for each of 51 discrete altitudes. The approximation of Figures 40 and 
41 are used to reduce execution time and cove storage. The minimum gap con- 
sidered, Gq, is t/cos 9 which is less than the smallest value of V versus 
inclination. When two latitudes (9^ and S^) are specified, the longitude 
viewed function V uses the minimum of that for the two extremes of latitude. 

When the inclination is such that two bands of longitude will be 
viewed, the second largest gap distribution is used (shown by a dotted line 
in Figure 37). This accounts for subdivision of the largest gap by the 
second "view" of the desired latitude. 



FIGURE 42. MAXIMUM GAP VERSUS ALTITUDE AND INCLINATION 



Maximum Gap 
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When a large range of altitudes is considered, there is a risk 
that none of the 51 discrete altitudes considered will fall in a forbidden 
(Inadequate coverage) region. The forbidden coverage (values of S) may fall 
between two discrete altitudes. This would result in a "clear** area on 
the screen with no warning to the user. If the space b tween the discrete 
altitudes is large enough for this to occur, evaluations are made for 
additional discrete altitudes between those 51 usually used. Any **forbidden" 
conditions encountered in the sub interval are plotted as if they occurred 
at the nearest usual discrete altitude. 

This coverage model has been exercised on a large number of cases 
and has proven to be an excellent means of indicating which orbits have the 
desired Earth coverage patterns. 


Radiation Environment 


The IGOS radiation model involves a straight-forward use of 
data published in Reference 33. These data are log-log plots of equiv- 
alent 1 Mev electron fluence per day /ersus altitude for several shielding 
densities at inclinations of 0, 30, 60, and 70 degrees. IGOS computes 
the maximum altitude versus inclination display by linear interpolation 
of the log-log data. 

Orbit-to-Sun Precession 


The orbit-to-Sun precession lines are plotted by solving Equation 
31 for the inclination as a function of altitude and the allowaMe precession 
rate 


I 


-1 

m cos 


K 

P 



(37) 



*f i 




9d 


where the allowable precession rate is computed from the allowable orblt/Sun 
precession angle, 0 , and time, T , (both user entered data), and the angular 

SO 60 

rate of the Earth's motion about the sun, U) , 

* es* 


(j) a U3 ± (J /T 

p es so so 


(38) 


Plotting Equations 37 and 38 results in two lines of altitude versus 

inclination, between which the orbit will process relative to the sun 

less than 0 in time, T . 

so so 


Orbit Decay 


The IGOS decay model operates by nvimerical (4th-order Range - 
Kutta) Integration of 


-1 



(39) 


where dl/dh is the derivative of orbit lifetime as computed in Reference 
34. Using the technique of Reference 34, for circular orbits, 


dL 

dh 


B • f (t,h) 
c s 


, dfe(h) 
dh 


. fjd) 


where 


B 

c 

f3(t.h) 

fj(h) 

fjd) 


the spacecraft ballistic coefficient, 
the solar activity factor, 
the normalized lifetime, and 
the inclination factor. 


(40) 
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IGOS contains the data for f^ and £^, aa shown In Figures 45 
and 46, which were taken directly from Reference 34. 

The solar activity factors requires estimation of future solar 
activity. Predictions to 1988 were obtained from MSFC and supersede the 
values given in Reference 34. The solar activity factor is computed from 

f (t.h) - f.(t/ (41) 

s t 

where the exponent X is 

X - .89027 (3 + 3.5 a - .5 , aud 

(Y - (h - 360)7240 

with h expressed in km. The time factor ^|.(^) plotted in Figure 47. 

IGOS Implements the functions shown in Figures 45, 46, and 47, 
by performing linear interpolation and the appropriate unit conversions. 

The differential equations. Equations 39 and 40, require the 
derivative of fj^(h). This is computed numerically from 

df^(h) fj(h + dh) - fj(h - dh) 

dh Tdii 


with dh set to 0.5 n. ml. Integration steps are computed by 


dt 


Min 


fat . 'I 

max I 

/ or ► 

"‘^'^max/i; 


(44) 


With dt « 0.2 years and dh “0.1 n. mi. Thus, no integration 
max max 

step <»vrceds either 0,2 years or 0,1 n, mi. Initial conditions are entered 
by the user and integration is stopped when time equals the specified mission 
duration, or altitude decays to the value represented by the botto^j of the 
IGOS plot. 
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FIGURE 47. SOLAR ACTIVITY TIME FACTOR f^(t) 
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Launch Vehicle Injection Error Model 

IGOS requires data for the dispersions in injection altitude 
(semi-major axis) as a function of orbit altitude and inclination. For all 
current conventional launch vehicles, the inclination dispersions are 
negligible for ordinary purposes and are not considered by IGOS. 

In its current version, IGOS contains only a Scout error model. 
The other vehicles utilize closed loop guidance and the resulting errors 
are small enough as to be, in general, negligible. If the injection errors 
of these vehicles are of interest for specific sensitive applications, they 
could be included by creating injection error tables similar to those in- 
cluded for the Scout. The data can be obtained from the launch vehicle 
project offices. 

Scout Error Model 


The Scout Project Office and LTV suggested Reference 35 as 
the best source of Scout accuracy data. This report includes altitude, 
velocity magnitude, and flight path angle dispersions. These data were 
combined t"o obtain the equivalent circular orbit semimajor axis dispersions 
shown in Figure 48. IGOS performs linear interpolation on the data to 
obtain the one sigma altitude dispersions as a function of nominal orbit 
altitude. 

Launch Vehicle Performance Model 

The IGOS launch vehicle performance model accesses a mass 
storage data file with a separate file record for each vehicle. To add 
a vehicle to the IGOS data file, the following data are required: 

(1) horizontal inertial velocity capability of the 

booster vehicle for an east-launch from ETR with a given 
pay’oad (Pj^) to a given altitude (A^). Up to 20 payloads 
and 5 altitudes may be used. 

(2) V4^, the delta-velocity of the fourth stage 

(Scout vehicles only) carrying payload, P^, and 

(3) the Indices of the sites from which the vehicle may 
be launched. 



Altitude Dispersion, n* tni 








I 
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These data may be obtained from published launch vehicle performance 
documents, or created using the NLVP- Performance program described in 
Reference 2 (or other similar programs). 

When the IGOS user specifies a particular launch vehicle, a 
set of vehicle performance tables are created from the performance data 
file. At this time, the range of altitude and inclination of interest 
have been defined. These tables consist of; 


(1) the payload capability for a circular orbit at 

altitudes and inclination, 1^, 

(2) the launch site used in computing ®nd 

(3) the launch azimuth used in computing 

The program computes 41 inclinations and 9 altitudes. The 9 
altitudes are distributed evenly across the range of the IGOS plot. For 
the inclinations, a set of 25 predetermined inclinations, shown in Table 5, 
are examined and those that lie in the range of the I'X)S plot are selected. 
The set of 41 Inclinations is then filled by interleaving additional incli- 
nations chosen evenly across the range of the IGOS plot. This technique is 
used to insure sufficient detail in the plots near the regions of plot 
discontinuities. Thus, the Inclinations have been chosen near the launch 
site latitudes and near the inclinations for which the preferred launch site 
changes with inclination. 


A set of launch site data is also stored. These data consist 
of launch site latitude and the launch azimuth constraints for each site. 

The payload table, is computed by computing the payload 

capability from each launch site permissible for the vehicle. If more 
than one site is permissible, the site giving the greatest payload is 
used. The direct ascent launch azimuth is computed. If this azimuth 
exceeds the launch azimuth constraint for that site, the constraint value 
is used and a plane change penalty is computed. 

With the payload table computed, the user then enters a 

specific payload value (PD) of interest. An altitude versus inclination 
contour is then plotted by linear interpolation at each inclination to 
compute the altitude corresponding to PD. 






CONCLUSIONS AND RECOMMENDATION'S 


The following summarizes Battelle's recommendations for further 
research based on the current status of the work for each task and the 
current t ends of NASA advanced studies activities. 

Task 1 


The computer codes, described In Table 2, have generally proven 
to be useful In the conduct of VSPRE studies. Several, however, have been 
widely up »d beyond the scope of VSPRE. Others could also, with further 
developTiu.' it , have similar continuing utility. 

EOMl'-I 


In addition to EOMP use on VSPRE at Battelle, and by NASA at 
MSFC, It has frequently been used by Battelle on both the DIGS (Delta Inertial 
Guidance System) project and NLVP (NASA Launch Vehicle Project). Some of 
these users also made several useful modifications to It: 

(1) Impact point dispersions (MSFC) 

(2) Metric output (MSFC) 

(3) Orbital launches for Shuttle upper stages (BCL/NI.VP). 

It Is recoiimended tnat EOMP- I be maintained and modified as necessary to 
maintain Its applicability In satisfying NASA booster dispersion analysis 
needs. 

IGOS 

This progran has been given wide visibility within NASA and 
has generated considerable Interest. 

The CO. Kept of IGOS, allowing quick response assessment of 
candidate orbi.es In defining mission requirements. Is of Interest to mission 
planners a well as launch vehicle and astrlonlcs specialists. The display 
of doma^ s of acceptable orbits actually Identifies candidate orbits which 
minimize the impact of injection dispersions. For defining astrlonlcs re- 
quirements In a ml’slon planning (orbit selection), the current IGOS Is only 
a preliminary cpvrational program. Several additional features have already 
been Identfled for inclusion In IGOS. These are discussed In the following 
pa cagraphs. 
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Specific Site Coverage^^# The IGOS Earth observation model presumes 
that the experimenter seeks frequent coverage of all longitudes in a given 
latitude band. For some missions, particular geographic locations are of 
interest; or, even with full longitude coverage, specific communication sites 
must be contacted for data down-links* An alternative coverage model could 
Indicate orbits which view a specific set of sites. 

Side Looking and Multiple Sensors. Many sensors view a portion 
of the Earth along a path displaced from the sub-satellite ground track. 

Also, many missions carry multiple sensors, each with a different set of 
viewing requirements. SEASAT, for example, carries multiple sensors, some 
of which are side looking. Straight-forward coding modifications to IGOS 
could enable it to handle .uich satellite missions. 

Elliptical Orblts^*^. The present version of IGOS considers only 
circular orbits. While extensionof the philosophy of the IGOS design volume 
to a three-dimensional space (apogee, perigee, and inclination) is straight 
forward, graphic presentation of three-dimensional spaces is difficult. Most 
likely, the best approach would be to continue to display a two-dimensional 
plot with the third-dimensional spaces is difficult. Two examples would be: 
apogee and perigee at a specified inclination, or eccentricity and inclination 
at a specified semi-major axis. Tf extensive application of IGOS is contem- 
plated, such modifications should be considered. 

Shuttle Performance Calculations^^. With the exception of the 
launch vehicle performance model, all IGOS features are applicable to the 
Shuttle as well as conventional launch vehicles. If it were assumed that 
the users mission would determine the orbital parameters for an entire 
Shuttle flight, the Shuttle (and its upper stages, if any) could be added 
to the vehicle data file. However, most missions will use only a portion of 
the Shuttle capability. Within this framework a number of types of contours 
could be drawn by IGOS to aid the user in selecting an orbit. For example: 


^These efforts have been included in a new study, at Battelle funded by MSFC. 



• The Shuttle user charges to any orbit, 

a The extra delta-V required to reach any orbit from 
the lowest cost or "nearest" currently scheduled 
Shuttle orbit placement, 

• Regions of common acceptability for experiment 
pallet sharing with other payloads. 

All of these require access by IGOS to data files which contain current 
Shuttle schedules and other Shuttle payload data. Many groups whithin NASA 
are considering particular aspects of the Shuttle scheduling and planning 
problem. Some of these are tending toward automated data management which 
will result in data files similar to those that would be needed for IGOS. 

It is recumtnended that the concept of IGOS receive continued 
funding. This will insure that the capability will exist to provide a 
quick espouse user oriented tool for orbit selection. With this tool, 
communications between the needs of the experimenter community and those 
responsible for NASA's transportation system capability could be conducted 
in a manner compatible with the needs associated with achieving the most 
efficient mission designs. 

Other Computer Codes 

The following computer programs have been written to support various 
VSFRE investigations. The source deck cards remain in the Battelle and/or 
MSFC files. They require minor modifications due to recent computer center 
operating systems changes if they are to be used in the future. 

EOMP- I I . This program was developed in 1972 to the point that, if 
adequate mission requirements are specified, it would be extremely useful 
on an operational basis. Since that time, effort has centered on the speci- 
fication of mission requirements, resulting in the development of IGOS. During 
this time, NASA's attention has concentrated on astrionics requirements for 
the Shuttle era. Application to Shuttle upper stage astrionics would require 
a major revision of EOMP-II. For this reason, it isrecommended that no 
further support be given EOMP-TI until interest justifies revisions to it. 

Monte Carlo . This program has served as a useful statistical 


analysis tool on a number of studies at BCL and MSFC. It is simple to use 



and requires no significant maintenance. It is recommended that the Monte 
Carlo code be maintained in its present form. This would require minimal 
effort. 


ABBACUS . ABBACUS was developed to aid in the evaluation of com- 
puter technology studies at MSFC. The funding for these technology studies 
is being revised. Continuation of the development and maintenance of 
ABBACUS should be commensurate with NASA's overall high reliability computer 
technology efforts. 

Spun Stage Simulation. This is a small program useful in the 
study of a number of problems related to spinning stages. In addition to 
Scout, for which it was developed, spun stages are used on Delta and Atlas/ 
Centaur and are now under consideration for a number of Shuttle era vehicle 
configurations. It is recommended that the spun stage program be maintained. 
This would require minimal funding. 

Task 2 


A number of astrionics requirements studies has been performed 
under Task 2. Most recently, these studies have been directed at the require- 
ments for high reliability computers. As mentioned in the discussion on 
ABBACUS recommendations, this particular study has been given low priority 
during the restructuring of NASA's combined high reliability computer funding. 

More generally, the evolution of the Shuttle upper stage program 
will create the need for astrionics-mlssion studies similar to those per- 
! srmed on VSPRE. It is recommended that the capability to perform these 
studi’..' be maintained. 


Task 3 

Task 3 was created when there was a high interest in the Advanced 
Small Launch Vehicle (ASLV), a major, evolutionary, upgrading of Scout. Since 
then Interest in ASLV has diminished. Concurrently, a number of less expensive 
modifications to Scout guidance have been proposed, analyzed by VSPRd and LTV, 
and implemented. It is anticipated that future improvement to the Scout may 
be considered. However, these will be of limited scope and could be analyzed 
under an effort similar to the current VSPRE Task 2. It is recommended that 
Task 3, be eliminated. 


With the creation of IGOS, a tool now exists to establish the 
accuracy requirements for expendable launches. It is recommended that 
an exercise be conducted to use IGOS to evaluate a set of baseline missions 
which will utilize NASA vehicles (conventional and Shuttle with upper 
stages). The evaluation should begin by considering the fundamental mission 
objectives such as coverage requirements or orbit lifetime. Using IGOS 
the domains of acceptable orbits can be established. The size of the domain 
may then be used for comparison to a set of launch system accuracies to 
establish which systems would be likely candidates and which would benefit 
from modification. These modifications could involve either accuracy 
improvement or cost reductions through using less accurate systems. If 
any system seems likely to benefit from modification, a cost study should be 
conducted to establish the cost/benefit when the system modification costs 
are included with per launch system hardware cost changes. 
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APPENDIX A 
IGOS PROGRAM ACCESS 

To use the Interactive Graphics OrbJt Selection program, it is first 
necessary to establish communications with the host computer. For the 

Battelle Computer System, conmunications can be obtained by dialing the 
Battelle data lines as follows: 

(1) (614) 421-7000 

(2) (614) 421-710^ 

(3) (614) 421-7040 

In the event difficulty is encountered, the status ^f the Battelle Computers 
may be checked by calling ('14) 291-9766 for a recorded message. 

Successful connection will be indicated by output of khe date and time and 
a request to "Login". In response to the login request type: 

LOGIN , VS PRE , ORBIT , SUP , N 

The terminal will indicate successful login by the request "COMMAND-". The 
following underlined commands should then be entered: 

COMMAND - ATTACH. PRG. COVER. ID°VSPRE 
COMMAND - EFL. 50000 
COMMAND - ETL.lOO 
COMMAND - PRG 

These commands connect the program file, extend the central memory field length 
to SOOOOg, extend the time limit to lOOg seconds, and begin execution of the 
program. 

Program execution begins with the display of a heading and an inquiry to 
determine if a graphics terminal is being used. 

The analyst is then able to exercise the program at his discretion using the 
command repertoire described later in this Appendix. 
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When execution of IGOS is completed (either by the IGOS 
command END, or by encountering a fatal execution error) the response 
•’COMMAND-" will be displayed on the bcreen. Reentering "PRG" will restart 
execution. If desired, a record of all user commands can be sent to the 
printer at Battelle (Building 13) by entering 

DISPOSE , XAPP99 , PR*IAH. 

If the session is complete enter LOGOUT to terminate the session. 

ent:aing-program commands 


A summary of all IGOS program command options are listed 
in Table A-1. Ihe commands are grouped into the following functional 
areas : 

e AXIS GENEBATION 
e COVERAGE DISPLAYS 
e VEHICLE PERFORMANCE 
e RADIATION DISPLAY 
e DECAY AND ORBITAL ACCURACY DISPLAY 
e SUN-SYNCHRONOUS DISPLAY 
e GRAPHICAL PLOT LABELING 
e TABULAR DISPLAYS 
e EXITTING PROGRAM 
e MISCELLANEOUS 

Whenever user input is required, the program issues a request by ringing the 
terminal bell and displaying two dashes (i.e., " - - " ). 

Associated with several commands are data variables that define 
the operation initiated by the command. The variable mnemonics and initialized 
values are given in columns 2 and 3 of Table A-1, respectively. Whenever a 
command request is executed, the current values of defined parameters are 
used to perform the operation. The parameter values can be changed in either 
of the following two ways 

(1) Individually - The parameter mnemonic is entered followed 
by the desired .rajue. For example, the maximum inclination 
for operation code "A" would be changed by entering - 


A-3 


MAX1NC*125.0« This modification capability can be used 
for the variables in Table A-1 marked with an 
(2) List • The parameters are Included with the command 
entry as a list. Parameter values are separated by 
commas and the relative position in the command string 
determine which data values are modified. For example, 
the maximum inclination variable for operation code "A" 
is the fourth (4th) parameter in the list. Therefore, 
the maximum inclination value would be changed by 
entering - A. ...125.0 . 

Whenever a data parameter is specified, the parameter's value remains 
equal to that value until again modified. All values can be reset to 
their initialized value by entering the RESET command. 

With a plot displayed on the screen, all commands are entered in 
a small rectangular box on the lower left of the screen. After the first 
command has been entered, it becomes difficult or impossible to read what 
has been typed. Therefore, the user must be attentive to his typing ability. 
If a:i error occurs and is detected in data entry, the user is informed by a 
multiple ring of the terminal bell. The command should then be reentered. 

No input entered by the user should cause abnormal program termination. 
However, if the error is not detected (l.e. , a wrong numerical value entered 
for some parameter) which results in an incorrect plot, there is no capability 
to erase the results. The terminal, that the IGOS program uses is a "non- 
refresh" screen type. Whatever appears on the screen is permanent until the 
entire screen is erased. (See Page A-7, Terminal Hardware.) If an error does 
occur, the user should request a new axis and all other data (i.e., coverage, 
radiation, launch vehicle payload curves, etc.) to make the plot complete. 

The user only need enter the basis commands since all data parameters (except 
where the error occurred) are set from the previous plot. 

The use of the above features are best illustrated in the example 
analysis session provided later in Appendix B, Sample Graphics Terminal 
Session. 



TABLE A-1. IGOS COMMAND OPTIONS 
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TEACH FEATURE 

The Interactive Graphics Orbit Selection program has a feature 
referred to as "TEACH" to help the participating analyst use the Interactive 
program. At any time In an exercise, all command options available at that 
current point are displayed merely by entering the "?" symbol. If questions 
then arise concerning the use of any specific command, Instructions for that 
command can be called for by entering the command followed by the "?" syndiol 
(l.e., COVER?). If a graphical display Is currently on the screen, a page 
request to erase the screen will be Issued. Permission Is given by entering 
any character and striking the return key. 

The teach feature Is Illustrated In Figures A-1 through A-3 which 
was taken directly from the terminal screen. 



COMPUTER HARDWARE 


The IGOS program was developed to use a Tektronix 4010/4012 series 
Interactive graphics terminal. The display screen Is a "non-refresh" type 
storage tube. The display area Is S.SO Inches by 6.25 Inches and consists 
of 1024 addressable points on the horizontal (X) axis and 781 points on the 
vertical (Y) axis. The terminal may operate In an alphanumeric (printing) 
mode or a graphic (plotting) mode. 

The keyboard Is a standard ASCII keyboard with several additions. 

Ihere are extra keys for erasing the screen and for causing hard copy 
generation. There are also twe thumbwheels for positioning the cross-hairs 
on the screen. Once positioned, striking any key on the keyboard will 
transmit the position of the Intersection of the cross-hairs to the computer. 

The cross-hairs are turned on and off by the user’s program. (l.e. , L C(M1AND). 

The hard copy unit Is used to generate a permanent (paper) copy of 
what Is displayed on the display screen. Approximately 5 seconds are required 
to generate a hard copy. The process does not affect what Is displayed on 
the screen. However, while the screen is sweeped during the copy process, no 
information can be displayed. If any Information comes from the computer during 
the sweep, It is lost. 
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(Continued) 
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Although the IGOS program utilizes the Tektronix terminal, 
it is possible to modify the program to operate on other terminals with 
similar features with a minimum of effort. This has been done to give 
demonstrations on Computek 400/15 terminals located at NASA Headquarters. 

The Battelle and NASA Headquarters terminals are currently 
driven by a Control Data CDC6400 computer system. A program version to 
use a Univac 1108 computer system has been developed and delivered to 
MSFC. The program is written in standard FORTRAN and should be adaptable 
to other computer systems with a minimum of conversion effort required. 
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APPENDIX B 


SAMPLE IGOS WORK SESSION 



APPENDIX B 


SAMPLE GRAPHICS TERMINAL SESSION 


The use of the Interactive Graphics Orbit Selection program is 
best illustrated by the following sample terminal session. All commands 
available to the analyst are used in this sample session. The displays 
'^hown were taken directly from the terminal screen. Table B-1 explains 
what each page illustrates and serves as an index to the work session. 



B-2 
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PAGE 


lABLE B>1. ORBITAL DESICN WORK SESSION GUISE 


DESCRIPTION 


B*4 Certain entries must he made to enter the program. Command EFL extends 

the field length, ETL e.xtends the time limit, TERM defines the number 
of characters per line, and ATTACH connects the computer program. The 
PEG convoand initiates execution. 

- B*5 The program is inquiring about the terminal being used. A "Y" response 

allows graphical output. 

B*6 An axis request is made with altitudes ranging from ISO to 300 n mi 

and inclinations ranging from 10.0 to 120.0 degrees. 

B-7 The resulting axis from the command request on page B»6. 

B>8 A coverage model using the standard defaults is requested. 

B*'9 The coverage model is hatched. 

B«10 Vehicle number 1 (i.e., SCOUT) is requested. 

B*ll Payload curves for 223.0, 300.0, and 350.0 pounds are requested and 

labeled using the "L" command. 

B*12 The coverage model from page B-3 is included with the payload curves 

shown on page B-11. 

B-13 A sun synchronous line is requested. 

B-14 The sun synchronous line is superimposed on the display created 

on page B-12 . 

B'lS The results of a radiation model request with .25 cm of aluminum 

shielding and a 2000 day exposure. 

B-16 The results of the radiation model are cross "hatched. 

B-17 The results of the radiation model, coverage model and payload curves 

are superimposed on a single axis display. 

B"18 A "ZOOM" is done for altitudes 165.0 to 200 nautical miles and 50.0 

to 75.0 degrees Inclination. Then payload curves for 225.0, 300.0 
and 350.0 pounds are superimposed. 'Hiis region appears to represent 
a potential orbital area to view 35.0 degree latitudes at least once 
every 30 days . 
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B-3 


TABLE B-1. ORBITAL DESIGN HORK SESSION GUIDE 
(Continued) 


PAGE 


DESCRIPTION 


B-19 The results of the decay model for a 240 n ml altitude, SO degree 

Inclination, a ballistic coefficient of 220, a 1977 launch data, and 
a mission duration of 3 years. 

B-20 Ihe results of the decay model for a 240 n ml altitude, 65 degree 

Inclination, a ballistic coefficient ot 220, a 1977 launch date, and 
a mission duration of 1 year. 


B-21 Ihe results of the coverage model from page B-8 and the decay model 

on page 6-20 superimposed. Some interesting conclusions can be 
drawn from this figure. First note that the Initial altitude distri- 
bution (ticks to the left) Indicate the possibility of inadequate 
coverage near + and -2^ . Noting the decay of +3cx and +10' it can be 
seen that an injection into the forbidden band would not decay out 
of the band for the entire mission. The lower forbidden band is 
much narrower and the decay rate much faster. Therefore, an injection 
in or above this band would result in a brief period of insufficient 
coverage. Attention should be paid to the fact that a -3o injection 
decays to 100 n.tni by May, 1977. These characteristics might warrent 
consideration of a slightly higher nominal altitude and lower 
inclination, for example 260 n.mi at 33°. 


B-22-23 

B-24-25 

B-26-27 

B-28-29 

B-30-31 

B-32 


Example of the payload vs. altitude and inclination data table. 

Example of the launch site vs. altitude and inclination data table. 
Example of the launch azimuth vs. altitude and inclination data table. 
Example of Che altitude vs. payload and iiiclination data Cable. 

Sensor table illustration. 

Vehicle library index. 


connANO- CFL.seeee 
COnnAND> ETL. 100 
COrinAHD- TERM. 85 

COIWAND- ATTACH. PRG. COWER, Iff- USPRE 
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